Digital Signal Processing 56 (2016) 79-92 











et 
ELSEVIER 


Contents lists available at ScienceDirect z unian 
igna 
Processing 


Digital Signal Processing 





www.elsevier.com/locate/dsp 





An evidence clustering DSmT approximate reasoning method for more 


than two sources 


® CrossMark 


Qiang Guo*, You He, Tao Jian *, Haipeng Wang, Shutao Xia 


Research Institute of Information Fusion, Naval Aeronautical and Astronautical University, Yantai, Shandong 264001, China 





ARTICLE INFO 


ABSTRACT 





Article history: 
Available online 27 May 2016 





Keywords: 

Information fusion 

Evidence theory 

Approximate reasoning 
Dezert-Smarandache Theory (DSmT) 
Proportional Conflict Redistribution 6 
(PCR6) 





Due to the huge computation complexity of Dezert-Smarandache Theory (DSmT), its applications 
especially for multi-source (more than two sources) complex fusion problems have been limited. To 
get high similar approximate reasoning results with Proportional Conflict Redistribution 6 (PCR6) rule 
in DSmT framework (DSmT + PCR6) and remain less computation complexity, an Evidence Clustering 
DSmT approximate reasoning method for more than two sources is proposed. Firstly, the focal elements 
of multi evidences are clustered to two sets by their mass assignments respectively. Secondly, the convex 
approximate fusion results are obtained by the new DSmT approximate formula for more than two 
sources. Thirdly, the final approximate fusion results by the method in this paper are obtained by the 
normalization step. Analysis of computation complexity show that the method in this paper cost much 
less computation complexity than DSmT + PCR6. The simulation experiments show that the method in 
this paper can get very similar approximate fusion results and need much less computing time than 
DSmT + PCR6, especially, when the numbers of sources and focal elements are large, the superiorities of 
the method are remarkable. 


© 2016 Elsevier Inc. All rights reserved. 





1. Introduction 


With the rapid development of science and technology, sensor 
types become diverse. To solve the more and more complex prac- 
tical problems, information fusion from different multi sources has 
been drawn great attention by scholars in recent years [1-12]. Due 
to the influence of noise, uncertain signal and information pro- 
cessing has become an important research direction in the field 
of information fusion. Belief function theory (also called evidence 
theory) has played a key role in uncertain and even conflict in- 
formation processing. As a traditional evidence theory, Dempster- 
Shafer theory (DST) [13,14] is a general applied information fusion 
method. However, DSmT, jointly proposed by Dezert and Smaran- 
dache [12], beyonds the exclusiveness limitation of DST and espe- 
cially in highly conflict information cases, it can obtain more accu- 
rate fusion results than DST. Recently, DSmT has many successful 
applications, such as, Map Reconstruction of Robot [15], Decision 
Making Support [16,17], Target Type Tracking [18], Image Process- 
ing [19], Sonar Imagery [20], Data Classification [21-25], Clustering 
[26-28], and so on. Besides, neutrosopic theory [29-31] proposed 
by Smarandache is a novel effective uncertain information process- 
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ing method. However, the main problem of DSmT is that when the 
number of sources and focal elements increases, the computation 
complexity of PCR5 or PCR6 in DSmT framework increases expo- 
nentially [32]. 

There were some important methods for reducing the computa- 
tion complexity of the combination algorithms in DSmT framework 
since this problem can be treated in different ways: 1) reducing 
the number of focal elements [33-36], 2) reducing the number 
of combined sources [37], 3) reducing both the number of focal 
elements and the number of combined sources [38]. Applied math- 
ematics has drawn attention by many scholars [39,40]. Particularly, 
the very recent Evidence Clustering DSmT Approximate Reasoning 
Method based on Convex Function Analysis [11] proposed by Guo, 
He, et al can get very similar approximate fusion results with PCR5 
in DSmT framework and cost little computation complexity. Nev- 
ertheless, the approximate method in [11] is only for two sources 
and it is not associative in the fusion of multiple (more than 2) 
sources of evidences. 

For reducing the huge computation complexity of Dezert- 
Smarandache Theory (DSmT) for multi-source (more than two 
sources) complex fusion problems and get more similar approx- 
imate fusion results with PCR6 in DSmT framework (DSmT + 
PCR6), a DSmT Approximate Reasoning Method for More than Two 
Sources is proposed in this paper. In Section 2, the basics knowl- 
edge on DST, DSmT + PCR6 and evidential distance theory are in- 
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troduced briefly. In Section 3, a new form of PCR6 formula is given 
and the mathematical analysis of the new form is conducted. Then 
a new DSmT approximate reasoning method for more than two 
sources is proposed. Finally, Analysis of computation complexity 
shows that the method in this paper costs much less computation 
complexity than DSmT + PCR6. In Section 4, the simulation exper- 
iments show that the method in this paper can get very similar 
approximate fusion results and need much less computing time 
than DSmT + PCR6, especially, when the number of sources and 
focal elements is large, the superiorities of the method are remark- 
able. In Section 5, the conclusions are given. 


2. Basic knowledge 
2.1. Dempster-Shafer theory (DST) 


A discernment frame based on the Shafer’s model is defined as 
© = {01, 02, --- , On} which contains n exclusive elements. The mass 
assignments of evidences defined over the power-set 2? is defined 
by 


m(Xi):2 10,11 Xic2®? (1) 


If m(X;) > 0, X; is called the focal element. m; denotes the mass 
assignments of the ith source of evidence. The Dempster’s combi- 
nation rule is given by [13,14] 


2 xinx;=x, izj M1 (Xi) - M2(Xj) 





X)= , WXCO 2 
mps(X) cK c (2) 
K= , mX) mX) (3) 
Xi, XjCO.iAj 
XiNXj=ø 


where K denotes the conflict beliefs of evidences. However, when 
the conflict beliefs of evidences are high, the fusion results of 
Dempster’s combination rule are usually very unreasonable. For 
this reason, many combination rules were developed, especially, 
Proportional Conflict Redistribution1-6 (PCR1-6) rules in DSmT 
framework have many advantages and successful applications [12]. 


2.2. Dezert-Smarandache theory (DSmT) 


The discernment frame of DSmT based on the hyper-power set 
abandons the exclusiveness limitation of DST. The hyper-power 
set denoted by D® admits the intersection of the elements. For 
example, if there are two elements in discernment framework 
© = {61,0}, the hyper-power set is D? = {Ø, 01, 02, 01 U62, 0162}. 
The mass assignments of evidences defined over the hyper-power 
set D® is defined by 


m(X;):D° — [0,1], Xi eD? (4) 


The Proportional Conflict Redistribution (PCR) rules [41,42] are 
the combination rules in DSmT framework. PCR rules have PCR1-6 
rules and the difference of them is that the proportional redis- 
tribution way of the conflict beliefs. PCR1-5 rules are applied 
for the combination of two sources and among these rules, PCR5 
is considered as the most precise mathematical method. PCR6 
rule is usually applied for more than two sources fusion prob- 
lems. 

PCR5 rule for 2 sources is introduced as follows [41,42] 


m1@2(Xi) = D mı(Y) -m2 (Z) (5) 


Y,ZeG® and Y,Z#ø 
YNZ=X; 





mpcrs (Xi) 
mıg2(Xi) : ; 
+ © reo and igi O + R] 
~ | Xe co and Xi 4G s 
0, Xi=ø 


where G? can been seen as the power set 2°, the hyper-power 
set D? and the super-power set SÊ, if discernment of the fusion 
problem satisfies the Shafer’s model, the hybrid DSm model, and 
the minimal refinement 0"! of © respectively and where all de- 
nominators are more than zero and the fraction is discarded when 
the denominator of it is zero [41,42]. 

This paper is mainly for more than two sources fusion. PCR6 
rule for more than 2 sources is introduced as follows 


M1@2@.-.@s(X) 


= > my (Y1) x m2(Y2) x -++ x ms(Ys) 


Y1, Y2, ,YsEG® and Y1,¥2,~- ,Ys#ø 
Y1NY2N---NAYs=X 


(7) 


MconflictTransfer (X) 
= 
> ; 2 
= 2Z4,Z2.,Zs-1€G? k=1(i1,i2, ,is)EP,2,-- ,5) 
ZiXX,i€{1,2,---,s—1} 
(Ni zDnx=ø 


f [= (X) x mj, (X) x +++ X Mi, (X) X Mi, (Z1) X ++ X Mi, 3] 


an 





[mi (X) + mig (X) +--+ +i, (X)] 








mi, (X) + Mi, (X) + +--+ mi, (X) +My, (Z1) +++ +My, (Zs—k) 
(8) 
mpcr6(X) = M192@..-@s(X) + MconflictTransfer(X), 
XeG? and XØ (9) 


where G® denotes the general power set which can be seen as 
the same as 2°, D® or the super-power set S? in different cases; 
and P(1,2,---,s) denotes the set of all permutations of the ele- 
ments. Equation (7) denotes that the combination products of the 
intersections of the mass assignments. Equation (8) denotes that 
the proportional redistribution of the conflict beliefs of mass as- 
signments. 
Assume that s = 2, PCR6 rule is given by 


m1@2(X) = D 


Y1,Y2€G? and Y1, Y2#ø 
YıNY2=X 


mı(Y1) x m2 (Y2) (10) 


mi, (X) x mea] 

mi (X) + mi (Z2) 
mi, (X) x Mi, oH 

Mi, (X) + mi, (Z1) 


MconflictTransfer(X) = Mj, (X) - | 


+ mj,(X) - | 


mpcre(X) = D mı(Y1) x m2 (Y2) 


Y1,Y2€G? and Y1,Y2#ø 
YıNY2=X 


mi (X) x mi, (Z2) 
mmea 

mi, (X) x mi, (Z1) 
o] 

X#D (12) 


+m, 00- 


+ mj, (X) | 


XeG? and 
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PCR has many advantages, such as it provides the appropriate re- 
distribution of conflicting beliefs, and it can produce reasonable 
combination result even in high conflicting cases. Nevertheless, the 
main disadvantage of PCR6 is that its computation complexity is 
too large and this limits its widely practical application. 


2.3. Evidential distance theory 


Evidential distance theory is a good way to measure the dis- 
similarity of multi evidences. Two well known dissimilarity mea- 
sure functions based on evidential distance theory are introduced 
briefly; first is Jousselme dissimilarity measure function based 
on Jousselme’s distance which takes into account the cardinal- 
ity of elements; second is Euclidean dissimilarity measure func- 
tion which has little computation complexity and fast convergence 
speed. 

1) Jousselme dissimilarity measure function Sim ;(m1, m2) [37] 
is defined based on the Jousselme’s evidential distance [43] as fol- 
lows 





1 
Sim (m1, m2) =1— Vm — m2)" Dim, — m2) (13) 


Re 


where D =[Djj] is a |G?| x |G®| positively definite matrix, and 
Dij =|XiN Xj|/|Xi UV Xj] with Xi, Xj € G®. 

The advantage of Jousselme dissimilarity is that it considers 
both the mass and the cardinality of bba’s. However, it makes no 
difference between the bba’s consisting of all single elements’ mass 
assignments and the bba’s consisting of non specific elements’ 
mass assignments |16]. 

2) Euclidean ESMS function Simg (m1, mz) [37] is defined based 
on the Euclidean evidential distance as follows 





|G?| 


Sime mma) =1- -5 Dim —m2(X;)]° (14) 


where |G®| is the cardinality of G®. 

The difference of Simg(m,,mz) and Sim)(m;,mz) has been dis- 
cussed in [37]. The detailed discussion is omitted in this paper. 
Please see [37] if necessary. It is proved in [37] that Simg(m1, m2) 
has many properties, including fastest convergence speed. So it is 
adopted in this paper as the dissimilarity measure of the approxi- 
mate method with DSmT + PCR6. 


3. An evidence clustering DSmT approximate reasoning method 
for more than two sources 


3.1. A new form of PCR6 formula for more than two sources 


It can be drawn from Equation (5)-(7) that the fusion proce- 
dures of the PCR6 rule have two main steps, first is the calculation 
of ™1@2@...@s(X) and second is the calculation of MconfictTransfer(X)- 
For reducing the computation complexity of our approximate rea- 
soning method, a new form of the PCR6 formula for more than 
two sources is given as follows 

1) The focal element of evidences from multi sources are de- 
noted by 


mı SOG: x}. eG? 
mz : {X?}, X? EG? 


ms: {Xf}, X$ € G? 


2) The first step fusion results of evidences are calculated by 


1 
MPCR-new (X) 








1 
=m, (X;,) 
3 ee x Mi, (X2) x +++ x me] 
s 1 ‘ 2 ; 7 
X2 X3 „XS EGO mi, (Xi) + Min (Xp) + + mi, (Xj.) 
Ig sige Mig 
if X}, NOXE NO- NX =O, X=Xi, 
else X =X}. NX} N NX? 
2 
Mbcr-new(X) 
2 
= mi, (X7,) 
y [ee x mi, (X?) Xe x a 
. 1 g 2 A 2 
do ys ee Min XP) mig (XB) + FM, OE) 
ty E as 


if X? NX) N---N Xi =O, X=X?, 
else X =X} NX} N Xf 


MPcR-new (X) 
= Mis (Xi) 


X1 X2 e, X671 eGo 
1 a Ts 4 





E (X$) x mi (X?) xX x o) 

mi, (X},) + mi, (XZ) H H mi KRO T 

if XP OXN NXT =H, X=X}, 

else X=X} NX] N NX] (16) 
3) The final fusion results are calculated by 


Mpcr-new(X) = Mocr-new(X) + Mocr-new(X) Laces Mcr-new(X) 
(17) 


3.2. Mathematical analysis of new form of PCR6 formula 


It can be drawn from Equation (16) that Micr-new(X). 
MŽer-new(X). “++, MÈcR-new(X) have the similar formula form. Due 
to the similar formula form, mber-new(X ) is analyzed first. 


1 
Mpcr-new(X) 
1 
=mMi, (Xi) 


x 


2 y3 e 
X? XÈ XR EG 


=mi (x1) x D [ms (XÈ) XX mi, (X;.) 





ee x mi, (X?) Xr x ar 
mi, (X$) + mi (XZ) +- + mi (X ) 





XÈ e, X$ eG 
3 s 
mi, (XÈ) 
x 2 T a 5 (18) 
acd mj, (Xi) + mi, (Xg) +--+ + mi, (X;) 
i2 
Let 

Mi, (x4) + Mi, (xÈ) +- + Mis (Xi) =a, 
mi, (X?) = x1, X2, Xn. (19) 
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Then 
aati 
X2 «G9 mi, (X$) + mj, (X?) +++ +i, (X) 
i2 
Z X1 X2 Xn 
~~ a+xı a+xX2 a+ Xn 
1 1 1 
=n—-ax + core (20) 
a+xı a+xX a+Xn 
Let 
fx) = , i=1,2,---,n, O<x% <1. (21) 
a+ Xj 


Because f (xj) is a convex function, the approximate convex func- 
tion formula [44] of )“"_, f (xi) is given by 











2 1 1 1 

YR) = = 

i a+x, da+x2 a+ Xn 
n 


(22) 





+A 
a+ (x1 +X2 +--+ +Xn)/Nn 


where A denotes the error of the approximate convex function 
formula, A > 0, A =0 iff x1 = X2 =- -- = Xn. Let xo = (X1 + X2 + 
---+Xn)/n, the convex function error analysis can be seen in [44] 
and the convex error function formula is given by 


= yy (xi — xo)? 
2(a + xo)? 


From Equation (23), the convex function errors are related to 
Yii — xo)? and ae If the average value Xo is large and 
X1, X2, ++ , Xn are concentrated, the errors are much smaller. 

Let n is the number of focal elements in the evidence, j is the 
source order number, j = 1,2,---,s and {xi}, i=1,2,---,n are 
the bba’s of the jth source’s evidence. 

The pseudo-code of Evidence Clustering method is given as Ta- 
ble 1. 


, X&<a+2xo (23) 


Table 1 
The pseudo-code of Evidence Clustering method. 





Input: 
The number of focal elements in the evidence: n 
The source order number: j= 1,2,--- ,s 


The bba's of the jth source’s evidence: gi), i=1,2,---,n 





1) For j=1,2,---,s 
{xj}, i=1,2,---,n are reordered to be {xi}, x} >x} > --- >x} in 
descending order. 
End 
2) For i=1,2,---,n 
Calculate the following function 


_ 15-0 ~The xp 


; (24) 
n—i 


F(x) 


If x! > fol), xi is forced to the set xi. 
Otherwise, xi and the following mass assignments are forced to the other 
set xi. 
End 

3) For j=1,2,---,s 
Calculate the sum and the number of xi. 
The sum of xÍ is denoted by xi and the number of xÍ is denoted by ki. 
Then the sum of A is denoted by x} =1- xi, the number of xi is 
denoted by ki =n- ki. 

End 





After the above Evidence Clustering method, the mass assign- 
ments are clustered to two sets, denoted by x} and Ri Besides, the 
sum and the number of two sets is calculated. 

Based on the Evidence Clustering method and the approximate 
convex function formula, Equation (16) can be transferred to 


1 
Mpcr-new(X) 


= Mj, (Xi, )? 


y {m:(X2) x +++ x mj, (X) 


XÈ XS eG? 
ig is 





e 
mi, (X$) + mi, (X?) +- mi (XÈ) 





X? G9 2 
i2 
xm; (X1 Y _ (x3 
wm (Xh) DO {mig (XR) x ++ x mi, (X}) 
XÈ s, XS eG? 
ig is 
Xa 
mi (Xj) + mis (XR) +- + mi, (XP) + X2 /kå 
X2 





b 
+ . x! 7 x3 ne í xs x2 z | 
mi, ( i) + Mi; ( p) bes + mi, ( i) + p/k 


3 
=mi (x1) ` 


X$ X$ EGO 
4 is 


{mi:,(x4) x +++ x mi, (X$) 





Xa 
1 3 S 2 j2 
%3 | (XE mis (X7, Hmi (XP XG ka ] 
x ) Mi, (Xi ) 1 z x2 3 


b 
3 
Xj, EG? mi (X; Hmi, (Xp )+=+mis (XP )+X5/kp | 


(25) 





Let 
mi, (Xi) + mig (Xj) +: + mis (Xi) + Xa/ka =a, 
mis (Xj) = x1 X2 Xie (26) 
Then 


3 x? 
D mi; (X;,)- pyi 7x3 ys 2 7p2 
mi, (X} ) + mi, (X}) ap “+++ mi, (XF ) + X46 /kG 


x X: Xi 
=x | ie a ee et a | (27) 
a+xı1 a+x a+Xn 








Based on the Evidence Clustering method and the approximate 
convex function formula, Equation (27) can be transferred to 


Xm, (X3): xa 
oh" ma (XE) + mig (XR) + mi XE) + X /ka 


x} EGO 
3 





N XG 
C ma (XÈ) + XG /ka + Xå /kå + mig (X3) +--+ mi, (XP) 
(28) 





Then Equation (16) can be transferred to 


1 
Mpcr-new(X) 


=Mi (x1)? > 


X$ o, XS EGO 
4 ts 


[m (X4) x ++ x mi (X) 
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x? x Xa 
mi, (X;,) + Xa/ka + Xa/ka + mig (Xf) +--+ + mi, (XP) 
3 
| X? x sa Xi 
: | mi, (X}.) + XB /kp + X2 /kå + mi, (X$) mi, (X$) 
| X? x Xa 
P mi (XE) + X3/k3 + X2/k2 + mj, (X$) mi, (X$) 
z x} 
b x 
mi, (X},) + Xp/kp + Xp /ke + mig (X) mi, (X) 























(29) 


Equation (29) can be transformed to the products of the sum of 
convex functions. So the mass assignments of focal elements are 
clustered by Evidence Clustering method for reducing the approxi- 
mate errors and the approximate formula of Equation (16) is given 
as follows 


1 
Mpcr-new(X) 
ees 
~ Mpcr-convex(X) 
1\2 
= mi, (X;,) 
J txt x- xt 
x 
mj, (X1) + t2/k2 + t3 /k3 +--+ t5/kS 
t?e[X2, X2}, k efk2,k2} aia aia A / 


t3e{X2, XP} ke e{k3 KP} 





te {XS, XS) kSe{kS KS} 
(30) 
where if ti = Xj, ki =k) and if t = XÍ, ki =k. 
Sag ; . j 
Similarly, the approximate formulas of MpcRr-new(X). 


j=1,2,---,s in Equation (16) are given as follows 


j 
Mpcr-convex (X) 





_m. (yi)? 
= mj;(X;,) 
| 
| Tigitets2,--,sy Mi (Xip | 
j l 

r jleg,2, st Mi Xi) + Digjet, 2, s] Mi Xip 
_m. (yi)? 
=m; ,(X;,) 


x 





Ix; ant 
D i #j,le[1,2,---,s] 
thetxt, Xi kek hy "i (Xi) + Lgjutett.2.- s E/K 
Kak!, if t!'=X!, 
Kom, if t!=X} 
(31) 


The Approximate Reasoning method is given as Table 2. 
3.3. Analysis of computation complexity 


Assume that s denotes the number of multi sources, and s > 2; 
n denotes the number of singleton focal elements in each evi- 
dence; there are only mass assignments of singleton focal ele- 
ments in multi source evidences, denoted by m(6;) > 0|6; € G? = 
{01,62,---,@n}; M, A, and D denote the computation complexity of 
multiplication, addition, and division for one time, separately. 

The computation complexity of PCR6 for more than two 
sources, denoted by opcre(n, s) is given by 


Opcre(N, S) 





=n-(s—1)-M+s-(n’—n)-[s-M+(s—1)-A+D]+n-A 


Table 2 
The pseudo-code of approximate reasoning method. 





Input: The number of focal elements in the evidence: n 
The source order number: j = 1, 2,--- ,s 


The bba’s of the jth source’s evidence: xi}, i=1,2,---,n 
the sum and the number of the two clustering sets of {x/}: X4, kå and Xj, ky, 
1) For j=1,2,---,s 

Calculate the following function 





j —m. (xi)? 
Mpcr-convex(X) = Mi; (Xj,) 
x% ay j lkl 
tle{X}, xb), ket kt} mi; (Xi) + Digjtett.2,--s) E/K 
Mak, if t!'=X} 
k'=k}, if t!=X} 


| 
Tigjtert,2,---,1¢ 





(32) 


End 
Calculate the convex approximate formula of Equation (17) as follows 


Mconvex(X) = MPcr-convex (X) + Mbcr-convex(X) +--+ Mbcr-convex(X) (33) 
Calculate the final approximate results by the normalization step as follows 


mconvex(X) 
y=. eae 34 
mea) Ž xece Mconvex(X) (34) 





=[n-(s—1)+s?-n- (n°! —1)]-M 
+ [s- (n° —n)-(s—1)+n]-A+s-(n°—n)-D (35) 


Proof. Assume that the multisource evidences are shown as fol- 
lows: 


mı(81) m4 (82) mı (6n) 
m2(@1) m2(02) m2 (6n) 

; (36) 
ms(61) ms(02) <-> ms(@n) 


The original PCR6 method can be divided into two parts as 
Equation (5), Equation (6) and Equation (7). 

1) First, the computation of Equation (5) can be represented by 
M1@2@...@s (0i) = Mı (8i) x M2 (8i) x --- x Ms (0i) (37) 


The computation of Equation (37) consists of (s — 1) times of mul- 
tiplication, denoted by (s — 1) - M. 

Because i = 1,2,---,n, the computation complexity of Equa- 
tion (5) of PCR6 method can be drawn as follows 


OPCR6-Equation (5)(1, $) =n- (s — 1) - M. (38) 


2) Then the computation of Equation (6) can be represented by 


MConflictTransfer (8i) 
Y1, Y2, ,Ys—1E€G? 
(QTI YONG =H 
mı (0i) x m2(Y1) x +++ x Ms(Ys—1) 
mı (6;) +m2(Y1) +----+ms(Ys—1) 


+ > my (6;) 
Y1,Y2,- ,Y¥s_1€G? 
NET Yone=0 
m2(6;) x m4(Y1) x --- x Ms(Ys—1) 
m2 (0i) +m (Y1) +--+ +ms(Ys_1) 


mı (0i) 
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ms (Gi) 


Jeete D 
Yi, Yoe ¥s-1€6? 
(Nii Y)NG=B 


ms(0i) x mı (Y1) X -++ X Ms_1(¥s—1) 








ms (0i) + mı (Y1) +--+ mMs—1(Ys—1) 
= X mæ 
Y1, Y2, ,Ys—1EG? 
mı (8i) x m2 (Y1) x +++ x Mms (Ys—1) 
mı (8i) 
mı (0i) + m2 (Y1) +--+ + mMs(Ys—1) 


mı (8i) x m2 (0i) x --- x Ms (0i) 
mı (6) + M2 (8i) +--+ + Ms (6i) 
+ > m2(6;) 

Y1,Y2,--- ,Ys_1€G® 

m (6) x m4(Y1) x +++ x Ms(Ys—1) 
m2(8i) + mı(Y1) +--+ +mMs(Ys—1) 
m2 (0i) x Mı (8i) x --- x Ms (0i) 

m2 (8i) + mı (0i) + --- + ms (0i) 








m2(0i)) 





+ 


i : 


Y1, Y2, .Y¥s-1€6® 


ms (0i) 


ms(0i) x mı (Y1) X <+- X Ms_1(¥s—1) 











ms (6;)) 
ms (0i) + mı (Y1) +-+: + mMs—1(Ys—1) 
ms (Gj) x mı (0i) x +-+ x Ms—1 (0i) (39) 
ms (0i) + m1 (8i) +--+ Ms_1 (Gj) 
One part of the Equation (39) is analyzed as follows 
mı (0i) x m2 (Y1) x --- x Ms(Ys—1) 
X mi) x 


Y1, Y2; ¥s_1€69 mı (0i) + m2 (Y1) +++: + Ms(Ys—1) 


sONGS mı (0i) x M2 (0i) x --- x Ms (0i) (40) 
mı (8i) + M2 (0i) +- - - + Ms (0i) 





The computation of Equation (40) consists of (n5™! — 1) -s times 
of multiplication, denoted by (n’~! — 1) - s - M; (nS! — 1) - (s — 1) 
times of addition, denoted by (n°7! — 1) - (s — 1) - A; and (n57! — 1) 
times of division, denoted by (n5™! — 1) - D. 

So the computation complexity of Equation (39) can repre- 
sented by s- (ns! — 1) - [s -M + (s — 1) - A + D]. 

Because i = 1,2,---,n, the computation complexity of Equa- 
tion (6) of PCR6 method can be drawn as follows 


OPCR6-Equation (6)(N, S) = S - (n° — n) -[s-M+(s—1)-A+D]. (41) 


3) Thirdly the computation of Equation (7) can be represented 
by 


mpcre (0i) = 192@...6s (8i) + MconflictTransfer (8i) (42) 


Because i = 1,2,---,n, the computation complexity of Equa- 
tion (7) of PCR6 method can be drawn as follows 


OPCR6-Equation (7) (1, 5) =n-A (43) 
Based on the above proof, the computation complexity of PCR6 
method can be drawn as follows 
opcre (N, S) 
= OPCR6-Equation (5) (1, S) + OPCR6-Equation (6)(N, S) 


+ OPCR6-Equation (7)(, S) 


=n- (s—1)-M+s- (n° -n)-[s-M+(s—1)-A+D]|+n-A 
=[n-(s—1)+s?-n-(n™!-—1)]-M 
+[s- (n° —n)-(s—1)+n]-A+s-(n°—n)-D (44) 














The computation complexity of the method in this paper, de- 
noted by ocy(n, s), is given as follows 
OGH(N, s) 
=s-n-25!.[s-M+(s—1)-A+s-D] 
+n-(s—1)-A+(n—1)-A+n-D 
=s’-n-25-!.M+[s-n-25"!-(s—1)+s-n—1]-A 
+ (s?-n-25' +n)-D (45) 


Proof. Assume that the multisource evidences are shown as fol- 
lows: 


m(@1) m4 (62) mı (0n) 
m2(891) m2(02) m2 (0n) 

. (46) 
ms(61) ms(02) -> ms (On) 


Each evidence is clustered by Evidence Clustering method and 
two sets of each evidence are generated. The sum and the number 
of two sets (sum, number) is obtained as follows 
(X4, ki) (Xi, kh) 

X2, k2) (X2, k? 
( a? a (X5 kp) (47) 


(XS, k) (X$, k$) 


The proposed method in this paper consists of Equation (32) to 
Equation (34). 

1) The computation of Equation (32) when j= 1 can be repre- 
sented by 


1 
™Mpcr-new (6) 
= m1 (6)? 


x 5 Px x- xt 
m4 (0i) + t2 /k2 + t3 /k3 + - - - tS /ks 
t2e{X2, X2}, K etk, k2} MEE / / 


t3e{X2, XP} Kee {k3 kp} 





te{ XS, XS}, kSe{kS, KS) 
(48) 
where if t = XÅ, ki =kd and if ti = xj, ki =k]. 


Because each of t?, t?,--- ,t° has 2 possibilities, 





txt x- xt 

my (0i) + t2/k2 + t3 /k? +--+ t8/ks 
t2e{ x2, x2}, K e{k2,k?} 1G) Fe} / / 
t3e{X2, XP}. Kee {k3,kP} 


t8e{XS, XEY, k efki KS} 


has 25—! items. 
Because i= 1,2,---,n, the computation complexity of Equa- 
tion (32) when j =1 can be drawn as follows 


OGH-Equation (32—j=1) (n, S) =N- E, [s-M+(s—1)-A+s-D] 
(49) 
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Table 3 
The numerical examples of computation complexity comparison. 
The increment of The method Numbers of division Numbers of addition Numbers of 
numbers of sources multiplications 
10 PCR6 method 1x10! 9x 10" 1x 10! 
The proposed method 10240 92250 102400 
50 PCR6 method 5 x 10°! 2.45 x 10” 2.5 x 10°3 
The proposed method 5.63 x 1016 2.76 x 1018 2.81 x 1018 
100 PCR6 method 1 x 10102 9.9 x 1018 1 x 10104 
The proposed method 1.27 x 10?2 1.26 x 1034 1.27 x 1034 
200 PCR6 method 2 x 10202 3.98 x 10704 4 x 10704 
The proposed method 3.21 x 10% 6.40 x 10% 6.43 x 1064 
2) Because Equation (32) when j= 1,2,---,s has symmetry, 1) The mass assignments of focal elements in each evidence are 


the computation complexity of Equation (32) can be drawn as fol- 
lows 


OGH-Equation (32)(N, $) = s-n- 2°"! -[s-M+(s—1)-A+s-D] (50) 
3) The computation of Equation (33) can be represented by 
OGH-Equation (33)(1, S) =n-(s—1)-A (51) 


4) The computation of Equation (34) can be represented by 








OGH-Equation (34)(N, S) = (n—1)-A+n-D (52) 
Based on the above proof, the computation complexity of the 
proposed method can be drawn as follows 
OcH(N, S) 
=s-n-25!.[s-M+(s—1)-A+s-D] 
+n-(s—1)-A+(n—1)-A+n-D 
=s*-n-25-!.M+[s-n-25!-(s—1)+s-n—1]-A 
+ (s?-n-25!+n)-D (53) 














It can be drawn from Equation (35) and Equation (45) that the 
computation complexity of PCR6 is almost proportion to s2-n’ and 
the computation complexity of the method in this paper is almost 
proportion to s*- 25. Analysis of computation complexity show that 
the method in this paper cost much less computation complexity 
than DSmT + PCR6, especially when the number of sources and 
the focal elements is large. 

Although the above theoretical analysis shows that the pro- 
posed method has higher computational efficiency than the origi- 
nal PCR6 method, the numerical examples are also given to verify 
this conclusion. Assume that there are only 10 singleton focal ele- 
ments which have mass assignments in hyper-power sets, denoted 
by G? = {61, 62, --- , 010}. The numbers of division, numbers of ad- 
dition and numbers of multiplications based on PCR6 method and 
the proposed method with the increment of numbers of sources 
are shown as Table 3. It can be drawn that the proposed method 
has more higher computational efficiency than the original PCR6 
method especially when the number of sources is large. 


4. Simulation experiments 
4.1. Simple cases of only singleton focal elements 


Example 1. If there are 3 evidence sources, assume that only 
singleton focal elements have mass assignments in hyper-power 
sets, denoted by Ge = {01,62,---,05}, k = 1,2,3. The mass as- 
signments in each evidence are m! = {0.1, 0.1, 0.3, 0.3, 0.2}, m = 
{0.2, 0.3, 0.05, 0.3, 0.15}, m? = {0.1, 0.05, 0.4, 0.35, 0.1}. 


clustered to two sets by the Evidence Clustering method, denoted 
by xq = {03,04,05}, xp = {01,02}; xq = (01,62, 04,05}, x5 = {03}; 
x3 = (04, 05}, XE = {01, 02, 03}. 
Then 

I el 1 1 = 
Xl =m} (03) + m4 (04) + m} (05) = 0.8 
X} =m! (0%) +m} (02) = 0.2 
X2 = mî (01) + m3 (62) + m4 (04) + m2 (05) = 0.95 
X? = m4 (03) = 0.05 
X2 = m3 (04) + mè (05) = 0.75 
XÈ = m? (61) + m3 (62) + m3 (03) = 0.25 (54) 





2) The convex function approximate fusion results are calcu- 
lated by Equation (32) 


mpcr-convex (1) 
_ mO- X-X 
m} (01) + Xa/3 + Xå /2 
my (01)? - X? - X? 
m! (01) + XP /2 + X} /3 
= 0.0169 


m} (01)? - X2. XP 
mj (01) + X¢/3 + X?/3 
mj (61)? - XÈ - X 
m1 (61) + X?/2+ X2 /2 











mpcr-convex (62) 
O m- X X 
~ m3 (62) + X3/3 + X3 /2 
my (02)? - X? - XÈ 
© mh (62) + X2/2 + X3/3 
= 0.0169 


m} (02)? - XZ - Xp 
m3 (62) + X¢/3 + X?/3 
m} (62) - XÈ - X 
m} (62) + Xp /2 + X2/2 








Mpcp-convex (3) 
m4(63)? - X2- X3 
~m3(03) + Xq/3 + Xq/2 
m3 (63)? - Xp - Xp 
m} (03) + X?/2 + XP/3 
= 0.1120 


m4 (63)? - X2 . XP 
m3 (03) + Xz/3 + X}/3 
m} (03)? - XÈ - X2 
m1 (63) + X2/2 + X32 /2 











Mbcr-convex (94) 
m1 (04)? - X2- X2 
~ m} (04) + X2/3 + X2 /2 
m1 (04)? - X? - XP 
m} (04) + X?/2 + X}/3 
=0.1120 


m} (04)? - X2 . X? 
m} (04) + Xa/3 + X}/3 
m} (04) - XÈ - X2 
m} (64) + X2/2 + X2/2 
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Table 4 
The fusion results and Euclidean similarities of the methods. 





1 2 


3 4 





The fusion results by 
our method 


0.0908, 0.1445, 0.1197, 0.4044, 0.2406 


The fusion results by 0.0909, 0.1451, 0.1194, 0.4045, 0.2401 0.2227, 0.1451, 0.1323, 0.2255, 0.2744 
PCR6 method 
Euclidean similarities 99.94% 99.92% 


0.2229, 0.1445, 0.1327, 0.2250, 0.2749 


0.2573, 0.2766, 0.1327, 0.2380, 0.0955 0.0778, 0.3109, 0.2648, 0.2380, 0.1085 
0.2570, 0.2768, 0.1323, 0.2384, 0.0954 0.0780, 0.3111, 0.2640, 0.2384, 0.1083 


99.95% 99.93% 





1 
Mpcr-convex 95) 
m! (05)? - X2. X2 
mi (05) + X2/3 + X2 /2 


mi (05)? - X2- XP 
mi (65) + X2/3 + X/3 
m} (05)? - X2 - X3 ms (0s)? - X? - Xa 
mi(0s) + X2/2+XP/3 mk (05) + XP/2+ X3/2 
= 0.0572 (55) 








Similarly, 


Mece-convex = [0.0572 0.1118 0.0047 0.1118 0.0348] 


56 
Mece-convex = [0.0183 0.0051 0.1875 0.1526 0.0183] 8) 


Then the convex approximate result is obtained as follows 
Mconvex = {0.0923, 0.1337, 0.3042, 0.3764, 0.1103} (57) 

3) The final approximate results are obtained by Equation (34) 
mou = {0.0908, 0.1315, 0.2991, 0.3701, 0.1085} (58) 

4) The fusion results of DSmT + PCR6 are obtained as follows 
MpsmT+Pcre = {0.0909, 0.1322, 0.2984, 0.3702, 0.1083} (59) 


The Euclidean similarity between mcy and mpcre is obtained by 
Equation (14) 


Ecu = 99.93% (60) 


From the above results of this example, the Euclidean similarity 
of the method in this paper can remain over 99.9% with DSmT + 
PCR6. It shows that the accuracy of the method in this paper is 
very high. 


Example 2. Assume that there are the same 3 evidence sources 
as Example 1. Ge = {01,02,--- ,@5}, k=1, 2,3 denotes the hyper- 
power sets of the 3 evidence sources. The mass assignments 
in 3 evidences are m! = {0.1,0.1,0.3,0.3,0.2}, m? = {0.2, 0.3, 
0.05, 0.3, 0.15}, mê = {0.1, 0.05, 0.4, 0.35, 0.1}. Move the mass as- 
signments of each focal element of m? one position backward at 
one time and 4 new evidences are obtained by 


m? = {0.1, 0.1, 0.05, 0.4, 0.35 
m’? = {0.35, 0.1, 0.1, 0.05, 0.4 
m? = {0.4, 0.35, 0.1, 0.1, 0.05 
m? = {0.05, 0.4, 0.35, 0.1, 0.1 


’ 


(61) 


’ 


koja ka a ore ee 


Each new evidence m? and the original evidence m! and m2 
are combined to obtain the fusion results by DSmT + PCR6 and 
the method in this paper. Then Euclidean similarities of two fu- 
sion results are obtained by Equation (14). The fusion results and 
Euclidean similarities of two methods are shown as Table 4. 

As shown in Table 1, the Euclidean similarities of the method in 
this paper with DSmT + PCR6 all remain over 99.9% and change 
little. It shows that the method in this paper has not only very 
high accuracy, but also high performance stability. 


4.2. Complex cases of both singleton focal elements and multiple focal 
elements in evidences 


Example 3. Assume that there are 3 evidence sources and the 
hyper-power sets are denoted by GE = {01 U 62, 01, 02, 03, 04, 05, 
8s, 87}, k = 1, 2,3. The mass assignments of evidences are denoted 
by m! = {0.2, 0.1, 0.1, 0.1, 0.3, 0.1, 0.05, 0.05}, m? = {0.05, 0.2, 
0.3, 0.05, 0.25, 0.01, 0.1, 0.04}, m? = {0.1, 0.1, 0.04, 0.4, 0.15, 0.1, 
0.1, 0.01}, the method process in this paper is given as follows 

1) The mass assignments of focal elements in each evidence is 
clustered to two sets by the Evidence Clustering method as fol- 
lows 


= {91 U 62,04}, x$ = (01, 02, 03, 05, 06, 07} 


1 
a 
x2 = {01 U 62, 01, 02, 03, O4, O6, 07}, xê = {05} (62) 


x2 = {03,04}, x} = {01 U 02, 01, 02, 05, 06, 07} 

Then 
Xx} =0.5, X}=0.5; Xł=0.99, X =0.01; (63) 
X2 =0.55; X?=0.45 


2) The convex function approximate fusion results are calcu- 
lated by Equation (32) 


Mice-convex = {0.0792, 0.0250, 0.0250, 0.0250, 0.1479, 
0.0250, 0.0072, 0.0072} 


macr-convex = {0.0071, 0.0771, 0.1442, 0.0071, 0.1093, 
0.0003, 0.0244, 0.0047} 


Mace-convex = {0-0257, 0.0257, 0.0049, 0.2297, 0.0510, 
0.0257, 0.0257, 0.0003} 


3) Then the normalized approximate convex results are ob- 
tained as follows 
Mconvex = {0.1014, 0.1157, 0.1576, 0.2370, 0.2791, 
0.0462, 0.0519, 0.0111} (65) 


4) For the singleton focal elements 63, 64,05, 06,67 which are 
not involved in the multiple focal elements, the approximate con- 
vex results of these focal elements are unchanged. For the multiple 
focal element 6; U 62 and the singleton focal elements 6; and 6, 
the convex approximate fusion results should be changed as fol- 
lows 

mMiekconvex ti U 62) has the wrong proportional mass assign- 
ments from the following items 


m! (6, U 02) @ m? (01) Q mM? (01), m! (01 U 62) Q m?(62) Q M? (62), 
m! (6, U 02) ® m? (01 U 02) ® M? (01), m! (01 U 02) D m2(6}) 
Q m3 (6; U62), 
m! (01 U 62) @ m? (01 U 62) & m? (02), m! (01 U 62) @ M (02) 
@ m? (01 U b2) 


which should be the mass assignments of mber-convex(61) and 


Mpcr-convex (2). Similarly, Mocr-convex (01 U 62). Mpcr-convex (1 U 
62) also have the wrong proportional mass assignment. So, the 
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Table 5 

Euclidean similarities of the results of the method in this paper with DSmT + PCR6. 
Experiment times 1 2 3 4 5 6 7 
Euclidean similarities 99.49% 99.46% 99.42% 99.41% 99.38% 99.46% 99.29% 





proportional mass assignments of 6; or 62 generated by 6; U 02 
are calculated as follows 


m-multiple! (6,, 62) = {0.0038, 0.0032} 


m-multiple” (01, 02) = {0.0003, 0.0003} (66) 
m-multiple? (01, 02) = {0.0015, 0.0018} 
Then 


mou (61 U 62) = mconvex(61 U62)— X` m-multiple! (61) 


j=1,2,3 
— > m-multiple/ (62) (67) 
j=1,2,3 
mou (91) =™Mconvex(@1) + >) m-multiple/ (01) (68) 
j=1,2,3 
mu (62) = mconvex(92) + J m-multiple/ (62) (69) 
j=1,2,3 


The final approximate results are obtained as follows 


Moy = {0.0906, 0.1214, 0.1628, 0.2370, 0.2791, 
0.0462, 0.0519, 0.0111} (70) 
4) The fusion results of DSmT + PCR6 are obtained as follows 


mpsmt+Pcr6 = {0.0858, 0.1232, 0.1678, 0.2380, 0.2789, 
0.0441, 0.0515, 0.0107} (71) 


The Euclidean similarity between mcy and mpcre is obtained by 
Equation (14) 


Ecu = 99.46% (72) 


It can be drawn in the above experiment results that in the 
complex cases of both singleton focal elements and multiple focal 
elements in evidences, the Euclidean similarities of the method in 
this paper with DSmT + PCR6 also remain over 99%. It shows that 
the method in this paper is also fit for the complex cases. 


Example 4. Assume that there are the same evidence sources as 
Example 3 and GP = {04 U6, 01, 62, 63, 04, 05, O6, 07}, k= 1, 2; 3 de- 
notes the hyper-power sets. Move the mass assignments of m? one 
position backward at one time to procedure 7 new evidences. 

The new 3 evidences are combined by DSmT + PCR6 and the 
method in this paper. Then Euclidean similarities of two fusion re- 
sults are obtained by Equation (9). Euclidean similarities of two 
methods are shown as Table 5. 

As shown in Table 2, in the complex cases of both singleton fo- 
cal elements and multiple focal elements in evidences, the method 
in this paper has not only very high accuracy which remains the 
Euclidean similarities with DSmT + PCR6 over 99%, but also very 
high performance stability. 


4.3. Associative analysis 


Example 5. Assume that there are 3 evidence sources and the 
hyper-power sets are denoted by Ge = {01,02,--- , 45}, k= 1, 2,3. 
The mass assignments in each evidence are the same as Example 1 
which are m! = {0.1, 0.1, 0.3, 0.3,0.2}, m? = {0.2, 0.3, 0.05, 0.3, 
0.15}, m? = {0.1, 0.05, 0.4, 0.35, 0.1}. 


1) The fusion sequence m! — m? —> m? is applied. The fusion 
results are from Example 1 as follows 


Mou = {0.0908, 0.1315, 0.2991, 0.3701, 0.1085} 


2) The fusion sequence m! —> m? —> m? is applied. The fusion 
process is given by 

(1) The mass assignments of focal elements in each evidence 
are clustered to two sets by the Evidence Clustering method, de- 
noted by 


1 = {63, 04,05}, x, = {01,62}; 

x3 = {04, 05}, xP = {01, 62, 03}; (73) 
2 = {01, 62, 04, 05}, x? = {03} 

Then 

Xxl =0.8, X}=0.2; X2=0.75, KP=0.25; (74) 
X2=0.95,  X?=0,05 


(2) The convex function approximate fusion results are calcu- 
lated by Equation (32) 
m-convex} (01) 

my (01)? - X2 - XZ 
m} (61) + X¢/2 + X2/3 
mi (01)? - XP - Xp 
m} (01) + Xp/3 + XP /2 
= 0.0169 


mi (01)? - Xp + X2 
m} (01) + XP /3 + X2/3 
my (01)? - X2 - X? 
m! (61) + X¢/2 + X?/2 











m-convex} (62) 
E m1 (02)? - X3 - X2 
m4 (62) + X3 /2 + X2/3 
m} (02)? - XP - Xp 
m1 (62) + X/3 + X?/2 
= 0.0169 


m} (02)° - XP - X2 
m1 (62) + X$ /3 + X2/3 
m} (02)? - X2 - X? 
m1 (62) + X2 /2 + X?/2 











m-convex} (03) 
_ 363)? - X XG 
m4(63) + X3/2 + X4/3 
m} (63)? - XP - X? 
m} (63) + Xp /3 + Xp /2 
= 0.1120 


m} (63)? - XP» XZ 
m3(63) + XÈ/3 + X2/3 
m} (63)? - X2 - X? 
m3 (63) + X2 /2 + X?/2 











m-convex, (44) 
O mO- XK? my (04)? - Xp -X2 
m4 (04) + X¢/2+Xa/3 m4 (04) + XP /3 + Xa/3 
m4(4) - Xp - X? m} (04) - X2 - X? 
© mi (4) + X2/3+X2/2  m}(04) + X2/2 + X?2/2 
= 0.1120 











m-convexs (05) 
mi (05)? - X2 - X2 
m1 (05) + X3 /2 + X7/3 


m} (05)° - Xp - X2 
mi (65) + XP /3 + X2/3 
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Table 6 
The mass assignments of the focal elements of four evidence sources. 
Highly conflict 81 62 03 04 
sources 
Sı x— 0.01 0.01 1—x-— 0.01 0.01 
S2 0.01 y— 0.01 0.01 1— y — 0.01 
S3 0.01 x— 0.01 1—x-— 0.01 0.01 
S4 y— 0.01 1— y— 0.01 0.01 0.01 





1.01 


o o 
© to 
œ% © 4 


Euclidean similarity 


o 
o 
= 





0.96 





Fig. 1. Euclidean similarities of the method in this paper with DSmT + PCR6. 


mi (05)? - Xp - XP 
ms(05) + XPp/3 + XP /2 


mi (05)? - X - X? 
mi (05) + X2 /2 + XP/2 





= 0.0572 (75) 
Similarly, 
m-convex? = [0.0572 0.1118 0.0047 0.1118 0.0348] 
m-convex? = [0.0183 0.0051 0.1875 0.1526 0.0183] ve) 
Then the convex approximate result is obtained as follows 
Mconvex = {0.0923, 0.1337, 0.3042, 0.3764, 0.1103} (77) 


3) The final approximate results are obtained by Equation (34) 
Moy = {0.0908, 0.1315, 0.2991, 0.3701, 0.1085} (78) 


By comparing the fusion process with Example 1, it can be 
drawn that the method in this paper has the associative property 
as the convex function approximate results of Equation (32) are 
unchanged with different fusion sequences. 


4.4. Cases of highly conflict evidence sources 


Example 6. Assume that there are four highly conflict evidence 
sources, and the hyper-power set of the four evidence sources are 
denoted by D? = {6;, 62, 63, 04}. The mass assignments of the focal 
elements of four evidence sources are shown as Table 6. 

Let x, y € [0.02, 0.98] and x, y increases from 0.02 to 0.98 by 
0.01 step. With the different values of x, y, Euclidean similarities of 
the method in this paper with DSmT + PCR6 are shown as Fig. 1. 
The average Euclidean similarity is also calculated as 99.17%. 

From Fig. 1 and the average Euclidean similarity, it can be 
drawn that under Cases of highly conflict evidence sources, the 
method in this paper can also get very high similar fusion results 
with DSmT + PCR6 which remains Euclidean similarities over 99% 
in the experiments and solve the highly conflict evidences fusion 
problems effectively. 


4.5. Convergence analysis 


Example 7. Assume that there are 4 evidence sources and the 
hyper-power sets of the evidence sources are denoted by G? = 
{01,62,---,@5}, k = 1,2,3,4. The mass assignments in evidences 
are m! = {0.01, 0.11, 0.35, 0.23, 0.15, 0.15}, m? = {0.01, 0.11, 0.23, 
0.35, 0.15, 0.15}, m? = {0.5, 0.3, 0.05, 0.05, 0.05, 0.05}, m* = {0.4, 
0.2,0.1,0.1,0.1,0.1}. Three methods are applied in this experi- 
ment to compare the convergence speeds. The methods are DSmT 
+ PCR6, the method in this paper and the DST combination 
method. 

First, fusion results of four evidences are obtained by differ- 
ent fusion methods. Then, the prior fusion results of the last time 
are fused with m?, m4 repeatedly. The fusion results by different 
methods for convergence analysis are calculated as Table 7. 

It can be drawn from Table 7 that 

1) The convergence speed of the method in this paper is faster 
than DSmT + PCR6 which converge to {0.7617, 0.1749, 0.0158, 
0.0158, 0.0158, 0.0158} in the 14th experiment and DSmT + PCR6 
converge to {0.7749, 0.1650, 0.0150, 0.0150, 0.0150, 0.0150} in the 
18th experiment. 

2) The convergence speed of the DST combination method 
is the fastest one among the three methods which converge to 
{1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000} in the 12th exper- 
iment, but this fusion result is strong converge to 6; and it is 
unchanged and wrong even if the following evidences focus on the 
other elements. Besides, the first two experiments’ results of DST 
are not reasonable for the high conflict evidences in these experi- 
ments. 

From the above convergence analysis, it shows that the method 
in this paper has good convergence speed and the convergence re- 
sults are also reasonable and effective. 


4.6. Monte Carlo simulation experiments in the case of random mass 
assignments in evidences 


Example 8. Assume that there are three evidence sources. The 
hyper-power sets of the sources are denoted by PpP = {01, 02, >, 
629} and only singleton focal elements have mass assignments. 
Monte Carlo simulation experiments are carried out 1000 times. 
The mass assignments of 20 focal elements in evidences are as- 
signed randomly each time. The fusion results of three random 
evidences are obtained by DSmT + PCR6 and the method in this 
paper. Then the Euclidean similarities of the method in this pa- 
per with DSmT + PCR6 are calculated and the computing time of 
methods are recorded as shown in Fig. 2, Fig. 3 and Table 8. (In this 
paper, all the simulation experiments are implemented by Matlab 
simulation in the hardware condition of Pentimu‘®) Dual-Core CPU 
E5300 2.6 GHz 2.59 GHz, memory 1.99 GB.) 

It can be drawn from Fig. 3, Fig. 4 and Table 8 that, in the 
case of random mass assignments in evidences, the method in this 
paper remains high Euclidean similarity and the average Euclidean 
similarity can reach 99.81%. Moreover, Euclidean similarity changes 
very little when the evidences are changed. Computing time of the 
method in this paper is much smaller than DSmT + PCR6. The 
simulation experiments results show that the method in this paper 
has strong practical meaning. 


4.7. Monte Carlo simulations in the case of increasing focal elements 
number and sources number 


Example 9. In this section, Monte Carlo simulation experiments in 
the case of increasing focal elements number and sources num- 
ber are performed. Assume that only singleton focal elements have 
mass assignments in hyper-power sets and the mass assignments 
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Table 7 
Fusion results by different methods for convergence analysis. 
Experiment DSmT + PCR6 The method in this paper Fusion results DST 
times 
1 0.5811, 0.2209, 0.0688, 0.0688, 0.0302, 0.0302 0.5796, 0.2218, 0.0681, 0.0681, 0.0312, 0.0312 0.0759, 0.8264, 0.0382, 0.0382, 0.0170, 0.0170 
2 0.6845, 0.2200.0.0272, 0.0272, 0.0205, 0.0205 0.6824, 0.2217, 0.0272, 0.0272, 0.0208, 0.0208 0.2326, 0.7599, 0.0029, 0.0029, 0.0008, 0.0008 
3 0.7284, 0.2022, 0.0177, 0.0177, 0.0171, 0.0171 0.7245, 0.2049, 0.0179, 0.0179, 0.0174, 0.0174 0.5048, 0.4948, 0.0002, 0.0002, 0.0000, 0.0000 
4 0.7494, 0.1870, 0.0159, 0.0159, 0.0159, 0.0159 0.7433, 0.1913, 0.0164, 0.0164, 0.0163, 0.0163 0.7728, 0.2272, 0.0000, 0.0000, 0.0000, 0.0000 
5 0.7607, 0.1775, 0.0154, 0.0154, 0.0154, 0.0154 0.7525, 0.1834, 0.0160, 0.0160, 0.0160, 0.0160 0.9189, 0.0811, 0.0000, 0.0000, 0.0000, 0.0000 
6 0.7670, 0.1719, 0.0153, 0.0153, 0.0153, 0.0153 0.7571, 0.1792, 0.0159, 0.0159, 0.0159, 0.0159 0.9742, 0.0258, 0.0000, 0.0000, 0.0000, 0.0000 
7 0.7705, 0.1688, 0.0152, 0.0152, 0.0152, 0.0152 0.7594, 0.1770, 0.0159, 0.0159, 0.0159, 0.0159 0.9921, 0.0079, 0.0000, 0.0000, 0.0000, 0.0000 
8 0.7725, 0.1671, 0.0151, 0.0151, 0.0151, 0.0151 0.7605, 0.1760, 0.0159, 0.0159, 0.0159, 0.0159 0.9976, 0.0024, 0.0000, 0.0000, 0.0000, 0.0000 
9 0.7735, 0.1662, 0.0151, 0.0151, 0.0151, 0.0151 0.7611, 0.1755, 0.0159, 0.0159, 0.0159, 0.0159 0.9993, 0.0007, 0.0000, 0.0000, 0.0000, 0.0000 
10 0.7741, 0.1656, 0.0151, 0.0151, 0.0151, 0.0151 0.7614, 0.1752, 0.0159, 0.0159, 0.0159, 0.0159 0.9998, 0.0002, 0.0000, 0.0000, 0.0000, 0.0000 
11 0.7745, 0.1653, 0.0151, 0.0151, 0.0151, 0.0151 0.7615, 0.1751, 0.0159, 0.0159, 0.0159, 0.0159 0.9999, 0.0001, 0.0000, 0.0000, 0.0000, 0.0000 
12 0.7746, 0.1652, 0.0150, 0.0150, 0.0150, 0.0150 0.7616, 0.1750, 0.0158, 0.0158, 0.0158, 0.0158 1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
13 0.7747, 0.1651, 0.0150, 0.0150, 0.0150, 0.0150 0.7616, 0.1750, 0.0158, 0.0158, 0.0158, 0.0158 1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
14 0.7748, 0.1650, 0.0150, 0.0150, 0.0150, 0.0150 0.7617, 0.1750, 0.0158, 0.0158, 0.0158, 0.0158 1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
15 0.7748, 0.1650, 0.0150, 0.0150, 0.0150, 0.0150 0.7617, 0.1749, 0.0158, 0.0158, 0.0158, 0.0158 1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
16 0.7748, 0.1650, 0.0150, 0.0150, 0.0150, 0.0150 0.7617, 0.1749, 0.0158, 0.0158, 0.0158, 0.0158 1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
17 0.7748, 0.1650, 0.0150, 0.0150, 0.0150, 0.0150 0.7617, 0.1749, 0.0158, 0.0158, 0.0158, 0.0158 1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
18 0.7749, 0.1650, 0.0150, 0.0150, 0.0150, 0.0150 0.7617, 0.1749, 0.0158, 0.0158, 0.0158, 0.0158 1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
19 0.7749, 0.1650, 0.0150, 0.0150, 0.0150, 0.0150 0.7617, 0.1749, 0.0158, 0.0158, 0.0158, 0.0158 1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
20 0.7749, 0.1650, 0.0150, 0.0150, 0.0150, 0.0150 0.7617, 0.1749, 0.0158, 0.0158, 0.0158, 0.0158 1.0000, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000 
0.127 ; > 
It can be drawn from the Table 3 that the computation complexity 
of PCR6 method is too large to make simulations when the num- 
ber of sources is big. So the maximum number of evidence sources 
> —©— Computing time of DSmT+PCR6 in cases 1-3 is set to be 5 and the number of focal elements is no 
2 ane ——— Computing time of the method in this paper more than 20 in the third case. In the fourth case, the number of 
= evidence sources increases with the constant number of focal el- 
D L . . 
£ 0.06 ements. Also for the large computation complexity of PCR6, the 
a maximum number of sources is considered as 10, and each evi- 
8 0.04- dence is composed of 5 focal elements. 
The simulation experiments are given as follows 
0.02; 1) Assume that there are 3 evidence sources, and the focal ele- 
ments of the hyper-power sets increase one focal element from 10 
0 i — i i 1 i 2 PE aera : : . Pi 
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Fig. 2. Computing time of the method in this paper and DSmT + PCR6. 
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Fig. 3. Euclidean similarity of the method in this paper with DSmT + PCR6. 


of the focal elements are assigned randomly each time. The exper- 
iments consist of four cases. In the cases 1-3, the number of focal 
elements increases with the constant number of evidence sources. 


Table 8 
Fusion results comparison in the case of random mass assignments. 


crement of the number of focal elements is recorded as shown in 
Figs. 4-5. 

2) Assume that there are 4 evidence sources, and the focal ele- 
ments of the hyper-power sets increase one focal element from 10 
to 50. The Euclidean similarities and computing time with the in- 
crement of the number of focal elements is recorded as shown in 
Figs. 6-7. 

3) Assume that there are 5 evidence sources, and the focal ele- 
ments of the hyper-power sets increase one focal element from 10 
to 20. For The Euclidean similarities and computing time with the 
increment of the number of focal elements is recorded as shown 
in Figs. 8-9. 

4) Assume that there are 5 focal elements, and the number of 
evidence sources increases one by one from 3 to 10. The Euclidean 
similarities and computing time with the increment of the number 
of sources is shown in Figs. 10-11. Especially, in order to explicitly 
present the increment of computing time, the ordinate scale is log 
instead of linear in Fig. 11. 

As shown in Figs. 4-11, 

1) The computing time of DSmT + PCR6 increases exponen- 
tially with the linear increment of the number of the focal ele- 
ments. However, the computing time of the method in this pa- 
per increases almost linearly with the linear increment of the 





Average Euclidean Max Euclidean 


Min Euclidean 


Average computing Max computing Min computing 





similarity similarity similarity time (s) time (s) time (s) 
DSmT + PCR5 0.1008 0.1153 0.1000 
The method in this paper 99.81% 99.98% 97.92% 0.0021 0.0306 0.0019 
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Fig. 4. Euclidean similarity of the method in this paper with DSmT + PCR6. 
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Fig. 5. Computing time of the method in this paper and DSmT + PCR6. 
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Fig. 6. Euclidean similarity of the method in this paper with DSmT + PCR6. 


number of the focal elements and remains very low computation 
time. 

2) The computing time of DSmT + PCR6 increases huge rapidly 
when the sources number increases. When there are 5 sources 
and 20 focal elements, the computing time of DSmT + PCR6 is 
123.2790 seconds. Whereas, in the case of 10 sources and 5 fo- 
cal elements, its computing time reaches 1213.9 seconds. It shows 
that when the number of evidence sources is large, DSmT + PCR6 
takes too big computation burden, and it is hard to satisfy the 
real-time requirement of intelligent systems. However, the com- 
puting time of the method in this paper increases little with the 
increment of evidence sources and it remains very little comput- 
ing time which can satisfy the real-time requirements of intelligent 
systems. 
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Fig. 7. Computing time of the method in this paper and DSmT + PCR6. 
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Fig. 8. Euclidean similarity of the method in this paper with DSmT + PCR6. 
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Fig. 9. Computing time of the method in this paper and DSmT + PCR6. 


3) In different conditions, the Euclidean similarities of the 
method in this paper with DSmT + PCR6 all remain over 99% 
and converge to very high value with the increment of the num- 
ber of the focal elements. It shows that the method in this 
paper can obtain very similar approximate fusion results with 
DSmT + PCR6, especially when the number of focal elements is 
large. 


5. Conclusions 


To get high similar approximate reasoning results with DSmT + 
PCR6 and remain less computation complexity, an Evidence Clus- 
tering DSmT approximate reasoning method for more than two 


Q Guo et al. / Digital Signal Processing 56 (2016) 79-92 91 


1.004 p 


—A— Euclidean similarity of the method in this paper 
1.002 - 


0.998 


0.996 


Euclidean similarity 


0.994 


0.992 








0.99 1 , , , , fl j 
3 4 5 6 7 8 9 10 


The number of sources 





Fig. 10. Euclidean similarity of the method in this paper with DSmT + PCR6. 


4 
Da Computing time of DSmT+PCR6 


Computing time of the method in this paper 




















= 
jo} 
N 
T 








Computing time (s) 
3, 





a 
9 
N 





10 i i i ; i i i 
3 4 5 6 7 8 9 10 


The number of sources 





Fig. 11. Computing time of the method in this paper and DSmT + PCR6. 


sources is proposed in this paper. Firstly, the focal elements of 
multi evidences are clustered to two sets by their mass assign- 
ments respectively. Secondly, the convex approximate fusion re- 
sults are obtained by the new DSmT approximate formula for more 
than two sources. Thirdly, the final approximate fusion results by 
the method in this paper are obtained by the normalization step. 
Analysis of computation complexity show that the method in this 
paper cost much less computation complexity than DSmT + PCR6. 
The simulation experiments show that the method in this paper 
can get very similar approximate fusion results and need much less 
computing time than DSmT + PCR6, especially, when the numbers 
of sources and focal elements are large, the superiorities of the 
method are remarkable. 

However, the mathematical analysis of convergence of the pro- 
posed method and the other evidence reasoning methods are still 
important tasks in our future research plans. We will also apply 
this new method in some real practical areas like radar signal com- 
binations in our future work. 
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